Parton Energy Loss in Two-Stream Plasma System 



Margaret E. Carrington 

Department of Physics, Brandon University, Brandon, Manitoba, Canada 

AND 

Katarzyna Deja 

National Centre for Nuclear Research, Warsaw, Poland 
AND 

Stanislaw Mrowczynski 

Institute of Physics, Jan Kochanowski University, Kielce, Poland 
and National Centre for Nuclear Research, Warsaw, Poland 

(Received January 6, 2012) 

The energy loss of a fast parton scattering elastically in a weakly cou- 
pled quark-gluon plasma is formulated as an initial value problem. The 
approach is designed to study an unstable plasma, but it also reproduces 
the well known result of energy loss in an equilibrium plasma. A two- 
stream system, which is unstable due to longitudinal chromoelectric modes, 
is discussed here some detail. In particular, a strong time and directional 
dependence of the energy loss is demonstrated. 

1. Introduction 

When a highly energetic parton travels through the quark-gluon plasma 
(QGP), it losses its energy due to elastic interactions with plasma con- 
stituents. This is the so-called collisional energy loss which for the equi- 
librium QGP is well understood, see the review p] and the handbook [2]. 
The quark-gluon plasma produced in relativistic heavy-ion collisions, how- 
ever, reaches a state of local equilibrium only after a short but finite time 
interval, and during this period the momentum distribution of plasma par- 
tons is anisotropic. Collisional energy loss has been computed for such a 

* Presented by K. Deja at the conference Strangeness in Quark Matter, Cracow, Poland, 
September 18-24, 2011. 
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plasma [3] but the fact that anisotropic QGP evolves fast in time, due to 
chromomagnetic unstable modes (for a review see [1]), has been ignored. 

We have developed an approach where the energy loss is found as the 
solution of an initial value problem. The approach is briefly presented in 
[3] where we also show that the formalism reproduces the well known result 
in case of an equilibrium plasma. In this paper we discuss in more detail 
a two-stream system, which is unstable due to longitudinal chromoelectric 
modes, and we demonstrate a strong time and directional dependence of 
the energy loss. 

Throughout the paper we use the natural system of units with c = h = 
ks = 1 and the signature of our metric tensor is (+, — , — , — ). 

2. Energy- loss formula 

We consider a classical parton which moves across a quark-gluon plasma. 
Its motion is described by the Wong equations [6] 



dx^(r) 

dr 
dp^(T) 

dr 
dQ a {r) 
dr 



«"(t), (1) 
gQ a (r)Fr{x(T))u u (T), (2) 
-9f abc n,(T)A^(x{T))Q c {T), (3) 



where r, x M (r), u^(t) and p M (r) are, respectively, the parton's proper time, 
its trajectory, four- velocity and four- momentum; F£ v and A% denote the 
chromodynamic field strength tensor and four-potential along the parton's 
trajectory and Q a is the classical color charge of the parton; g is the coupling 
constant and a s = g ,2 /47r is assumed to be small. We also assume that the 
potential vanishes along the parton's trajectory i.e. our gauge condition 
is Uh(t) Aa (x(t)) = 0. Then, according to Eq. the classical parton's 
charge Q c (r) is a constant of motion. 

The energy loss is given directly by Eq. ^ with {i = 0. Using the time 
t = jt instead of the proper time r and replacing the strength tensor Fa 
by the chromoelectric E a (i,r) and chromomagnetic B a (t,r) fields, Eq. ^ 
gives 

d ^l=gQ-B a (t,r(t)). V , (4) 

where v is the parton's velocity. We consider a parton which is very ener- 
getic, and therefore v is assumed to be constant and v 2 = 1. 

Since we deal with an initial value problem, we apply to the field and 
current not the usual Fourier transformation but the so-called one-sided 
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Fourier transformation denned as 

f(uj,k) = J dtj dW^^/frr), 



/(f, r ) = r m ±L f " e -'M-^)/( W) k), 



(5) 
(6) 



where the real parameter a > is chosen is such a way that the integral over 
uj is taken along a straight line in the complex uj— plane, parallel to the real 
axis, above all singularities of /(w,k). Introducing the current generated 
by the parton j a (i, r) = gQ a v5^ (r — vt) and using Eq. Eq. Q can be 
rewritten: 



r°™> du r d*k {uj _ Q)t 

L^IW 6 E a (.,k).v, 



dt 



(7) 



where u> = k • v. 

The next step is to compute the chromoelectric field E a . Applying the 
one-sided Fourier transformation to the linearized Yang-Mills equations, we 
get the chromoelectric field given as 

k) = -i^fiuj, k) [ U ji( u , k) + ei kl k k B l 0a (k) - uDi a (k)] , (8) 

where Bq and Dq are the initial values of the chromomagnetic field and the 
chromoelectric induction D l a {uj,\s) = £ ij (uj, \&)Ea{uj, k). The chromodielec- 
tric tensor e^(o;,k) equals 



e ij (uj,k) =5 ij + 



d 3 p 



df( P ) 



2oj J (2vr) 3 uj - k • v + i0+ dp k 



k- v 



UJ 



5 kj 



+ 



k k vi 



UJ 



where /(p) is the momentum distribution of plasma constituents. The color 
indices a, b are dropped because e(u, k) is a unit matrix in color space. The 
matrix £ l -?'(u;,k) from Eq. (IsJ) is defined 



E# (w, k) = -k 2 ^' + fe*jy + uj 2 e ij (uj, k). (9) 
Substituting the expression ^ into Eq. ([7]) , we obtain the formula 



dE(t) 
dt 



gQ a v 



duo r d 3 k 

-oo +l a 2^ J 

iuigQ a vi 



w — UJ 



+ WB 0a (k)-^ a (k) 
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When the plasma is stable, the poles of k) are located in the lower 

half-plane of complex u) and the corresponding contributions to the energy 
loss from these poles decay exponentially in time. The only stationary 
contribution to the energy loss is given by the pole uj = uj = k • v from 
the first term in square brackets in Eq. (10). Therefore, the terms which 
depend on the initial values of the fields, can be neglected, and Eq. (10) 
reproduces the known result for the energy-loss of a highly energetic parton 
in an equilibrium quark-gluon plasma [5]. 

When the plasma is unstable, the matrix S _1 (o;,k) contains poles in 
the upper half-plane of complex uj, and the contributions to the energy loss 
from these poles grow exponentially in time. In this case, the terms in 
Eq. ( 10 ) which depend on the initial values of the fields D and B cannot be 
neglected. Using the linearized Yang-Mills equations, the initial values Bq 
and Do are expressed through the current and we obtain 



* -^y^ww e ^'™'- k) (ll> 

x - (k j k k -k 2 ^)(S- 1 ) H (w,k)+ W we^(u,k)(S- 1 ) fc; (w,k)," 

-UJ — UJ 

where Cr comes from the averaging over color states of the test parton, 
Cr = 4/3 for a quark and Cr = 3 for a gluon. 



3. Two-stream system 



In order to calculate the energy loss, one must invert the matrix T, lJ (uj, k) 
defined by Eq. ([9]) and substitute the resulting expression into Eq. (11). 
For a general anisotropic system this is a tedious calculation. In the case 
of the two-stream system, which has unstable longitudinal electric modes, 
the chromodynamic field is dominated after a sufficiently long time by 
the longitudinal chromoelectric component. Consequently, we assume that 
B(w,k) = and E(w,k) = k(k • E(w,k))/k 2 , in which case the tensor 
E lJ (u;,k) is trivially inverted as 



(w,k) 



1 



k l k 3 



e L (w,k) = e«(w,k) 



uj 2 el(uj,^) k 2 
and Eq. (fTTj) simplifies to 

"°° +WT duo r 

2W (2vr)3 w 2 eL(u;5k) k2 



k l W 
"k 2 " 



dE(t) 
dt 



9 2 C R 



d 3 k 



-i(ui—ui)t jjj2 



0J 



— r + - 

— UJ UJ 



(12) 



(13) 



Eq. (13) gives a non-zero energy loss in the vacuum limit when el — > 1. 



Therefore, we subtract from the formula (13) the vacuum contribution, or 
equivalently we replace 1/el by 1/el — 1- 
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The next step is to calculate el(u, k). With the distribution function of 
the two-stream system in the form /(p) = (27r) 3 n <5( 3 )(p — q) + 5^ 3 ^(p + q 
where n is the effective parton density in a single stream, one finds [7J 



(u - w+(k)) [u + w+(k)) (uj - w-(k)) (w + w_(k)) 



(k-u) 2 )' 



(14) 



where u = q/-Eq is the stream velocity, and // 2 = g 2 n/2E Ci is a parameter 
analogous to the Debye mass squared. There are four roots to the dispersion 
relation £l(cj, k) = which read 



4(k) 



i 



k 2 (k-u) 2 +/i 2 (k 2 -(k-u) 



(15) 



± fiJ (k 2 - (k • u) 2 ) (4k 2 (k • u) 2 + ^ 2 (k 2 - (k • u) 2 )) 



It is easy to see that 0<w + (k)6l for any k. For k 2 (k-u) 2 > 2/i 2 (k 2 - (k- 
u) 2 ) the minus mode is also stable, < w_(k) G M.. However, for k • u ^ 
andk 2 (k-u) 2 < 2^ 2 (k 2 — (k-u) 2 ) one finds w 2 (k) < Oanda;_(k) imaginary. 
This is the well-known two-stream electric instability. 

Strictly speaking, the stream velocity u given by the distribution func- 
tion equals the speed of light. However, the distribution function should be 
treated as an idealization of a two-bump distribution with bumps of finite 
width. This means that the momenta of all partons are not exactly parallel 



or antiparallel, and the velocity u which enters Eqs. (14 15) obeys u < 1. 



Equations ( 13 ) and ( 14 ) determine the energy loss of a parton in the two- 



stream system. The integral over to can be computed analytically as it is de- 
termined by the six poles of the integrand located at oj = ±w+(k), ±o;_(k), 
Q and 0. The remaining integral over k must be done numerically. We use 
cylindrical coordinates (k = (kT,4>,ki,)) with the axis z along the streams. 
The parameters are chosen to be: g = |v| = 1, |u| = 0.9, Cr = 3. As 
discussed in [5] , the integral over k is divergent and it has been taken over a 
domain such that — /c max < < fc mffl and < hp < &max with /j max = 20//. 

In Fig. [T] we show the parton energy loss per unit length as a function 
of time for different orientations of the parton's velocity v with respect to 
the stream velocity u. The energy loss oscillates and manifests a strong 
directional dependence. 



4. Conclusions 



We have developed a formalism where the energy loss of a fast parton in 
a plasma medium is found as the solution of an initial value problem. The 
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Fig. 1. The parton energy loss per unit length as a function of time for several 
angles 9 between the parton's velocity v and stream velocity u. 

formalism determines the energy loss in an unstable plasma which contains 
modes that exponentially grow in time. The two-stream system has been 
studied in some detail. The energy loss per unit length is not constant, 
as in an equilibrium plasma, but it exhibits strong time and directional 
dependences. 

Acknowledgments 

This work was partially supported by Polish Ministry of Science and 
Higher Education under grants N N202 204638 and 667/N-CERN/2010/0. 

REFERENCES 
[1] S. Peigne and A. V. Smilga, Phys. Usp. 52, 659 (2009). 

[2] M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, 
2000). 

[3] P. Romatschke and M. Strickland, Phys. Rev. D 71, 125008 (2005). 

[4] St. Mrowczyhski, Acta Phys. Polon. B 37, 427 (2006). 

[5] M. E. Carrington, K. Deja and St. Mrowczyriski, |arXiv:1110.4846| [hep-ph]. 

[6] S. K. Wong, Nuovo Cim. A 65, 689 (1970). 

[7] St. Mrowczyhski, Phys. Rev. D 77, 105022 (2008). 



